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It has been proposed recently to observe the change in cosmological redshift of distant galaxies 
or quasars with the next generation of large telescope and ultra-stable spectrographs (the so-called 
Sandage-Loeb test). Here we investigate the possibility of observing the change in peculiar velocity in 
nearby clusters and galaxies. This "peculiar acceleration" could help reconstructing the gravitational 
\ potential without assuming virialization. We show that the expected effect is of the same order of 

, magnitude of the cosmological velocity shift. Finally, we discuss how to convert the theoretical 

f^*) ■ predictions into quantities directly related to observations. 
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INTRODUCTION 



One of the most outstanding issues in astrophysics and cosmology is the measurement of the total mass of clustered 
structures. The problem became especially crucial after it was realized that most matter in the universe is not directly 
visible through astronomical observations. Since direct investigation is impossible, a number of indirect techniques 
for mass determination have been developed over the years. The most frequently used methods rely on kinematic 
Qh measurements, where the velocity dispersion of some suitable class of test "particles" is used to infer the virial mass 
of the object under investigation. This class of methods were successfully applied, for example, to the rotation curves 
of spiral galaxies (which led to the first evidence of the existence of large amounts of dark matter) and to the velocity 
dispersion of elliptical galaxies. The largest clustered structures in the universe where the virial theorem is still 
applicable are the clusters of galaxies. The measurement of the velocity of galaxies in clusters led to further evidence 
in favour of dark matter, and so it did the observation of bremsstrahlung x-ray emission from free electrons of the 
. intracluster gas. 

Until now, dynamical methods (i.e. measurements of acceleration) have not been applied to the measurement of 
mass, essentially because of the extraordinary difficulties in detecting tiny velocity variations for astrophysical objects 
expected over time scales much smaller than the cosmic time. However, in the near future we will reach a level of 
technological development that might allow the observation of the dynamical effect of mass. In fact, measuring tiny 
drifts of the order of 1 cm/sec in the spectra of astrophysical objects over time scales of a few decades will be within 
reach of high-resolution, ultra-stable spectrographs, coupled to extremely large telescopes (with apertures of 40-60 
m) currently under planning fij. Possible applications of this kind of measurements have been proposed, for example, 
to detect extra-solar planets [2j], or for cosmological applications [3|, |4j, as a way to detect variations in the Hubble 
rate of expansion of the universe, thus investigating such issues as dark energy or the variation of physical constants 

sag. 

In this paper, we explore the possibility of using the time variation of the peculiar velocity of astrophysical objects 
(an effect that we refer to as peculiar acceleration [4]) as a way to dynamically determine the mass profile of clustered 
structures. This method does not assume virialization, although for simplicity here we assume a spherically symmetric 
gravitational potential. We find numerical estimates of the effect expected for some typical objects, under idealized 
conditions, and show that the predicted signal is of the same order of magnitude of the cosmological one. The 
determination of the peculiar acceleration is interesting also because it will act as a systematic noise source in the 
observation of the cosmological velocity shift. Since the detailed specifics for future experiments are still largely to be 
finalized, we leave a detailed comparison of the theoretical predictions with the observational settings to future work. 

II. PECULIAR ACCELERATION 

The peculiar acceleration of a particle at the spherical-coordinate position (r, a, 7) in a system centered on a cluster 



is 



a = V<J> = $ r f -\ r — $ a a + 7 7 (1) 

r sin 7 ' r ' 
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where the tilded quantities are unit versors and $ is the gravitational potential. We will always assume that Earth's 
local acceleration has been properly subtracted. The acceleration along the line of sight versor s is then 

a s = s • V$, (2) 

Let us assume now spherical simmetry for the potential, namely $ = <f>(r). Therefore the line of sight acceleration is 
as follows: 

a s = cosfl'^ r , (3) 

where 9' is as in Fig. [1] In the same figure we define R c as the cluster distance from observer, r as the particle distance 
from cluster's center and 9 the viewing angle. 

For small viewing angles 9 (i.e. for r <C R c ) we have that the peculiar acceleration along the line of sight is 

a s = cos6>'$ r w sin/3$ r |r=,R c 0/co S /3 (4) 
and assuming spherical symmetry we obtain 

n GM{r) 

a s pa sin/3 — — , (5) 

where r = R c 9/ cos/3. Finally, assuming a density profile p(r), we can write a s in terms of the viewing angle 9,R C , (3 
and the density profile parameters. 

Let us assume now for definiteness a Navarro- Frenk- White profile (NFW) [1, [§] for the dark matter halo responsible 
for the cluster gravitational potential 

p( r ) = j_(i P +j_y > ( 6 ) 

where r s = r v /c sets the transition scale from r -3 to r" 1 , c is a dimensionless parameter called the concentration 
parameter, p cr — 3Hq/8ttG is the critical density at the redshift of the halo (H being the Hubble constant), r v is the 
virial radius inside which the mass density equals A c p cr and S c is the characteristic overdensity for the halo given by 

5 C = (7) 

where 

C=[log(l + c)- 



l + c J 

In addition, A c is the nonlinear density contrast for a virialized object and enters the expression for the mass 
M v — 4/37rr 3 A c /9 cr . Its value depends on the cosmological model and assuming a ACDM we can set A c = 102 [To| . 
Then the mass associated with the radius r is 

r 

M(r)=M,C'(log(l + ^)- T ^ T ) (8) 

and consequently 

GM ( log(l+^) 1 \ , Q . 

Considering a time interval AT the velocity shift of the particle test due to the peculiar acceleration along the line of 
sight is s = Av — a s AT and turns out to be 

S = — CATSm P{ (Z)2 ^ (1 + ^) J- ( 10 ) 

Finally we obtain the velocity shift s(/3) 

n ,,cm . AT Af , r s /log(l + £) 1 v 
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with r — R C Q I cos j3. For a typical cluster value of C ~ 1, it turns out therefore that the typical shift for a galaxy 
cluster is of the order of 1 cm/sec , similar to the cosmological value at z » 1 [f| . The maximal value {(3 — > tt/2, 
r — > 0) is : 

„ „cmAT Mj, . r s 2 C 

So = - 44 ^i(^W 2" < 12 > 

In the following we estimate the strength of this velocity shift effect for different clustered objects, i.e. at different 
distances and scale. We will always assume AT = lOyr and let /? vary in (0, 7r/2). For /3 G (— 7r/2,0) the effect is 
clearly opposite in sign. 




Figure 1: Definition of geometric quantities. 



A. Virgo cluster 

The Virgo cluster is the nearest cluster of galaxies. It is also the most remote cosmic object with a physical 
connection to our own small group of galaxies, the Local Group. Its distance is about 15Mpc [ll|, corresponding 
to a redshift of 0.003 (assuming that Ho — 72km/s/Mpc): this ensures that the cosmological velocity shift signal 
does not come into conflict with the one due to peculiar acceleration (the cosmological velocity shift being of the 
order of 10 _2 cm/sec at this redshift). The virial radius is taken to be r v — 2.2Mpc corresponding to a virial mass 
of 1.2 x 10 15 M Q ; the concentration parameter is c = 4, leading to r s — 0.55 Mpc. In Fig. [2] the velocity shift as a 
function of (3 is shown for this cluster. In particular the curve is plotted for four different values of the viewing angle 
6. 
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Figure 2: Predictions for Av in cm/sec for Virgo (R c = 15Mpc, r v = 2.2Mpc , c = 4) for four values of 9 up to 9 max = 0.15 
equally spaced, starting from 8 m in = Omax/i (top to bottom). This value of 9 m in corresponds to a radius of 0.55Mpc. 



B. Coma cluster 

Coma is the most studied and best known cluster of galaxies, having three very appealing properties [HI, [12]: it is 
almost perfectly spherically symmetric, very rich (M v — 1.2 x 1O 15 M ), and close to our Local Group (R c = lOOMpc, 
z = 0.02 ; the cosmological velocity shift is of the order of 10 _1 cm/sec at this redshift). Assuming, as above, a NFW 
density distribution, the other parameters are: r v = 2.7 Mpc, c = 9.4 and r s — 0.29 Mpc [l2|. The velocity shift 
predicted for Coma is shown in Fig. [3l 

Coma 
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Figure 3: Predictions for Av in cm/sec for Coma (R c = lOOMpc, r v = 2.7Mpc , c = 9.4) for four values of 9 up to 9 max = 0.027 
equally spaced, starting from 6 min = 9 max /A (top to bottom). This value of 9 min corresponds to a radius of 0.675Mpc. 



C. Andromeda 

It may be interesting to apply this method on a typical galaxy scale. We chose to investigate as a representative 
case the Andromeda galaxy, located nearly at the center of our Local Group, at a distance of 0.8 Mpc from us. In 
order to fit the NFW profile to a galaxy, its parameters must be chosen in a way that does not directly reflect the 
immediate physical meaning. For Andromeda, we chose R s = 16kpc, r v = 200 kpc, and c = 12.5 (Q). Its mass is 
estimated to be roughly M = 10 12 M Q . Obviously in this case the test particles will be represented by stars or gas 
clouds, rather than galaxies. 
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Figure 4: Predictions for An in cm/sec for Andromeda (R c = 0.8Mpc, r v = 200kpc , c = 12.5) for four values of 6 up to 
dmax = 0.25 equally spaced, starting from m i„ = 6 max /4: (top to bottom). This value of 8 m in corresponds to a radius of 
7.5kpc. The signal increases for smaller radius. 



III. OBSERVING THE PECULIAR ACCELERATION 

Since observations are performed along the line of sight, it is practically impossible to obtain information on the 
galaxy position within the cluster, i.e. on the (3 angle. It is then useful to identify some observable quantities that 
can be directly compared to theoretical predictions. We propose two such quantities in the following subsections. 

A. Maximum velocity shift 

The predicted velocity shift reaches a maximum value along the line sight, Av max , that depends on the halo profile. 
By estimating Av ma x from a sample of galaxies along each line of sight, one can therefore constrain the cluster 
gravitational potential. Fig. [5] shows the behaviour of A 

Umax 3j function of 8 for Virgo, Coma and Andromeda. We 
note that each curve can be actually observed only for angles up to max . We also plot in Fig. [6] and Fig. [7] the same 
quantity Av max as a function of mass for a cluster and a galaxy, respectively. Clearly the signal is independent of the 
object distance for a maximal viewing angle 6 max ~ I /Re- 
in Fig. [8] and Fig. [9] we plot contour curves corresponding to the same Av max for 8 = 9 max /4 and different 
combinations of the concentration parameter c and mass. In principle, knowing the mass of the object from other 
observations, a measurement of the velocity shift can give indications on the concentration parameter. In particular, 
recently several authors have derived fitting formulae for the concentration parameter as a function of mass from 
TV-body simulations ([3, [H, HE])- We show a curve from one of these formulae in Fig. [HI Then, for clusters, a 
measurement of Av max at different mass scales would provide a test of the concentration parameter fitting formula. 
Clearly, instead of the maximum one could use the variance of the peculiar acceleration or another non-vanishing 
average. In the next subsection we use instead the full distribution to relate the theoretical prediction to observable 
quantities. 

B. Density distribution of galaxies 

Knowing the radial distribution of galaxies p g (r) (i.e. deriving it from the projected number density) one can 
compare the observed distribution N(As) (number of galaxies in the velocity shift bin As) with the prediction. 
The expected number of galaxies in an annulus AO and thickness A(i? c #tan/3) = R c 8(\/ cos 2 (3) A/3 around the 
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Figure 5: The maximum of the velocity shift due to peculiar acceleration as a function of the viewing angle 6 for Virgo (straight 
line), Coma (short-dashed line) and Andromeda (long-dashed line). The round points mark the value of 9 max lying on the 
correspondent curves. 
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Figure 6: The velocity shift signal as a function of mass for a Coma-like cluster, located at the same distance as Coma, for 
three different values of 8 separated by 8 max /3, starting from 8m ax /3 (top to bottom). 



cluster center is 



n(R c , 9, s)A0As = 2vp g ( ^] ^,A9A(3 (13) 
\cos p / cos^ p 



and therefore we expect in the interval As a number of galaxies proportional to 

RcQ \ R 3 J 2 Ae A.s 

cos [3 J cos 2 (3 (ds/d/3) 



n(R c , 9, s)A9As = 2ir Pg ( ) -£^0 ,^,^ (14) 



where (3 — /3(s) can be obtained by inverting ifTTj) . The histogram N(R C , 9, s) = nAs is therefore given by 

N(R c As) = 2n Pg (WL\ ~4^f As ^ 
\cosp(s) / cos / p(s) as 

Now Eq. lfT5|) contains only observable quantities. In principle, confronting this expected number N with the 
observed numbers one can reconstruct p(r). An example of N(s) is given in Fig. ifTTj) . 
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Figure 7: The velocity shift signal as a function of mass for an Andromeda-like galaxy, located at the same distance as 
Andromeda, for three different values of 9 separated by 9 max /3, starting from 9 max /3 (top to bottom). 
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Figure 8: Contour plot for the velocity shift of a cluster at 9 max /4: as a function of the concentration c and the mass M. 
Higher values correspond to darker regions and the contours levels are (left to right) Av — (0.1, 0.5, 2, 4, 6, 8, 10, 15)cm/s. The 
short-dashed line is the fitting formula by 



C. Comparing with the cosmological signal 

As already mentioned in the previous section, the cosmological signal is competitive with the velocity shift due to 
the peculiar acceleration. In Fig. [lOlthe cosmological signal for a ACDM model with fl\ = 0.7 is plotted together with 
the maximum of the velocity shift of a cluster; the redshift of their "equivalence" is very small (z ~ 0.02) depending 
on the mass of the cluster. In addition, there is a second redshift, around z ~ 2, at which the cosmological signal is 
negligible. However observations at such a distance could not be reliable, not only because of the difficulty of precise 
measurements at this redshift, but also because the zero-crossing of the function is cosmology-dependent. In principle 
of course the cosmological signal could also be subtracted from the peculiar one after averaging over a consistent 
number of galaxies. 
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Figure 9: Contour plot for the velocity shift of a galaxy at 8 ma x/4: as a function of the concentration c and the mass M. Higher 
values correspond to darker regions and the contours levels are (left to right) Av — (0.1, 0.5, 2, 4, 6, 8, 10, 15)cm/s. 
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Figure 10: The cosmological velocity shift as a function of redshift for a ACDM model (solid line) and the maximum of the 
velocity shift due to the peculiar acceleration for 9 = # maa; /4 and three different value of the mass of a rich cluster: M = 10 14 M© 
(short-dashed line), M = 5 • 1O 14 M (long-dashed line), M = 10 15 M© (long-short-dashed line). 



IV. CONCLUSIONS 



In this paper, we performed a first exploration of the magnitude of peculiar acceleration expected in clustered 
structures in the universe, over time scales of one or a few decades. Our study was motivated by the interesting 
prospects of observing velocity shifts of cosmic objects of order of a few cm/sec in a decade offered by forthcoming 
extremely large telescopes coupled to high- resolution ultra-stable spectrographs. Peculiar acceleration would be an 
interesting target for such future measurements, since, as we have shown, could give an independent measurement of 
the mass profile of clusters of galaxies, and possibly even of individual galaxies. We have also shown that measurements 
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Figure 11: Frequency histogram N(s) (expected frequency of galaxies in velocity shift bins) for Coma at 6 — 9 ma x/^ ■ Here we 
assumed that the radial distribution of galaxies is proportional to the halo profile. 

of the velocity shift for different mass scales could test the relation between the concentration parameter of clusters 
and the virial mass. 

Our results show that the expected signals for the objects we considered (the Coma and Virgo clusters, and the 
Andromeda galaxy) are comparable to the cosmological change in redshift, due to the expansion dynamics, expected 
for distant QSO's [1|]. Due to the high stability of Ly-a lines in QSO's, it has been recently argued that the cosmological 
effect could be observed in a time scale of about three decades [3, @| , giving important information on the composition 
of the universe, and in particular on the mysterious dark energy which drives its current accelerated expansion. It 
seems plausible that an effect of similar order of magnitude, predicted by our study, could then be also observable, 
although the objects used as test particles are different in nature. Our study is also useful to the extent of assessing 
the peculiar acceleration as a source of errors in the measurement of the cosmological velocity shift Q • 

Of course, we realize that the experimental and observational details of future instruments are still too blurry to 
make any clear statement on the actual feasibility of using peculiar acceleration as an astrophysical and cosmological 
tool. Many issues are still to be discussed before forecasting the results of actual observations, from the stability of 
the galactic spectra to the effect of local deviations from spherical symmetry. However, we believe it is exciting to 
entertain the possibility that peculiar acceleration might actually be observable in the not too distant future, and 
hope that more detailed investigations of its applicability in astrophysics might be stimulated by our study. 
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